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The phase diagrams survival-extinction for the Penna model with parameters: (mutations rate)- 
(birth rate), (mutation rate)- (harmful mutations threshold), (harmful mutation threshold)-(minimal 
reproduction age) are presented. The extinction phase may be caused by either mutational meltdown 
or overpopulation. When the Verhulst factor is responsible for removing only newly born babies and 
does not act on adults the overpopulation is avoided and only genetic factors may lead to species 
extinction. 
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I. INTRODUCTION 

A species becomes extinct when the last existing 
member of that species dies. Extinction therefore be- 
comes an extreme event after which no surviving speci- 
mens are able to reproduce and create a new generation. 
There have been at least five mass extinctions in the his- 
tory of life in which many species have disappeared in a 
relatively short period of geological time. 

The classical "BigFive" mass extinctions identified by 
Raup and Sepkoski Q were approximately 444 (End Or- 
dovician), 367 (Late Devonian), 255 (End Permian), 200 
(End Triassic) and 65 (End Cretaceous) million years 
ago. These mass extinctions may have been caused by 
impact events, climate changes, volcanism, 7-ray burst 
and plate tectonics or the mixtures all of them. Dur- 
ing End Permian extinctions — the worst in the Earth 
history — 96% of marine species and 70% of land ones 
(including plants, insects and vertebrates animals) were 
killed. The most recent, the most famous and most rapid 
mass extinction (65 million years ago, at the end of the 
Cretaceous period) wiped out all non-avian dinosaurs 
and it was caused probably by a cosmic collision with 
an asteroid or comet of ten kilometers in diameter [ij . 

We arc currently in the early stages of a human-caused 
mass extinction, known as the Holocene extinction event. 
Probably, up to twenty percent of all living species could 
become extinct within thirty years (by 2028). Wilson 
estimated that if current rates of human destruction 
of the biosphere continue, half of all species of life on 
Earth will be extinct in hundred years [ij. 

On the other hand even without external tragedies 
the fate of populations may be quite unpleasant when 
an organism's number exceeds the carrying capacity of 
its ecological niche. Such situation is termed as over- 
population. Thomas Malthus argued that if left unre- 
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stricted, human populations would continue to grow un- 
til they would become too large to be supported by the 
food grown on available agricultural land Q. He pro- 
posed that, while resources tend to grow linearly, popu- 
lation grows exponentially. At that point, the population 
would be restrained through mass famine and starvation. 
Malthus argued for population control, through "moral 
restraint" , to avoid this happening. As the population 
of a species exceeds the amount of available resources, it 
decreases, sometimes sharply since the lack of resources 
causes mortality to increase 

Could a population become extinct due to genetic fac- 
tors? The answer is probably positive. The Darwinian 
evolution theory [f| predicts genetic modification of in- 
dividuals genome via mutation process and slow change 
one species to another. This change (called speciation 
IE IE Hi) an d the newly created species will be accepted 
by Nature through the natural selection procedure. In 
this way mutations — usually regarded as harmful - 
may become crucial for new species creation and evolu- 
tion. However, sometimes too many genetic modifica- 
tions yield genetic death of the individual. Sometimes it 
happens for all individuals in a given population — we 
call that the mutational meltdown 91. 



II. THE PENNA MODEL 

The Penna model of biological ageing is devoted 
to reproduce single-species population dynamics for ge- 
netically heterogeneous individuals represented only by 
their genotype — the Nut long binary computer word. 
The population fate is assumed to be governed by accu- 
mulation theory [TlT |. which claims that random heredi- 
tary deleterious mutations accumulate over several gen- 
erations in the species genomes. These bad mutations are 
represented as "1" in the genome. The time is measured 
by discrete variable t. In each time step (which may be 
treated as a "season" and may have different meanings 
for different species), for each individual of age a, the 
number of bits set to one in the first a positions in the 
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genome is calculated. These "1" in position 1 < i < a 
are treated as active mutations. If the number of active 
mutations is greater or equal to the threshold value T an 
individual dies. 

Individuals compete among themselves for resources 
(for example food and territory) : each of them may be re- 
moved from the population with probability N(t)/N max , 
where N max represents the maximal environmental ca- 
pacity and N is the current number of individuals. This 
process is independent on individual age and health. The 
term N(t)/N max is termed a Verhulst's factor 01 The 
Verhulst factor in natural way limits exponential explo- 
sion of the population size to the infinity, to avoid the 
Malthus catastrophe. 

The random removing of organisms due to their com- 
peting for resources is not very well justified biologically: 
we feel that the well-fitted individuals should be more 
resistant for random removing — which mimics competi- 
tion — that ill- fitted ones. Thus, to keep the population 
finite and avoid accidentally killing the best-fitted organ- 
isms, the Verhulst factor may act only on newly born 
babies and not on adult ones [l^. Here, however, the 
idea of massacre of the innocents — proposed by Herod 
the Great Q — is applied to one season old children 
and not to children who are two years old and under. 

The population reproduces asexually employing 
cloning mechanism. If the individual is older than a min- 
imum reproduction age R but before loosing its fertility 
— which happens at age E — it clones itself with proba- 
bility b producing B offspring. The cloning process is not 
perfect: during replication parent's genome is exposed to 
harmful mutations ("0"— which occur with proba- 
bility m at M randomly selected positions in genome. 
The length of genome Nut restricts maximal age of indi- 
viduals. 

The extinction process within the frame of the Pcnna 
model was discussed rather seldom. In the literature 
known to us only the papers by Maksymowicz et al (lij . 
Bernardes 0, Pal H3 and Fehsenfeld et al [3 deal with 
species extinction. Some possible suggestions about mu- 
tational meltdown arc also available in Refs. [SIM EH- 
On the other hand the rapid survival rate decreases for 
semelparous organisms (like Pacific salmon j^]) or the 
influence of an overfishing on northern cod fish popula- 
tion |23[ were one of the first applications of the Penna 
model j2J|. 

In this paper we would like to draw a phase diagrams 
in the space of the Penna model control parame- 
ters which for initially large populations separate survival 
from extinction. 



III. THE RESULTS 

We start our simulation with No = 3 • 10 s individu- 
als with perfect genomes (N^u zeros) and with environ- 
mental capacity N max = 10 6 . We allow individuals to 
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(b) N max =10 6 , N =3 10 5 , T=3, R=8, E=32, B=5, m=b=1. 
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FIG. 1: A few examples of population size dynamics for dif- 
ferent mutation and birth rates (M, B) and probabilities (m, 
b). The subsequent pairs of curves have identical set of pa- 
rameters except of mutation rates which are M and M + 1. 
The Verhulst's factor applies either to (a) all individuals or 
(b) only to newly born babies. 



reproduce until their death (E = Nut = 32). 

In Fig. examples of the population density 

N(t)/N max are presented. Subsequent pairs of curves 
correspond to sets of (m, b) parameters (0.25, 0.25), (0.5, 
0.5) and (1, 1) from top to bottom, respectively The rest 
of parameters is fixed except the mutation rate M, For 
fixed set of parameters (T, R, E, B. m, b) the time evolu- 
tion of population density differs qualitatively for muta- 
tion rate M (survival) and M+l (extinction). Such and 
similar differences allow to construct phase diagrams pre- 
sented in Figs. EUSandH In Fig. □» the Verhulst fac- 
tor reduces only a birth probability b — > b[l — N(t) /N max ] 
and not the total number of individuals (see Fig. n^a)). 

In Fig. |21 the phase diagrams survival-extinction for 
various values mutation and birth rates (M, B) and prob- 
abilities (B, b) and for average numbers of mutations 
Mm and births Bb per "season" and per individual are 
presented. The phase transition border for average val- 
ues of number of mutations Mm and average number 
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14 
12 
10 
8 
6 
4 
2 




14 



12 



10 



(b) N max =10 6 , N =3 10 5 , T=3, R=8, E=32 
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(c) N max =10 6 , N =3 10 5 , T=3, R=8, E=32, m=b=1. 
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FIG. 2: The phase diagram survival-extinction for (a) various 
values of the mutation and birth rates (M, B) and probabil- 
ities (B, b) and (b) for average numbers of mutations Mm 
and births Bb per "season" and per individual. Below arc- 
like lines the species survival is possible. At and above them 
individuals die due to too many mutations pumped into their 
genome. On the right side of the vertical line the extinction 
process is caused by overpopulation, (c) The latter may be 
avoided when Verhulst's factor acts only on newly born ba- 
bies. 




FIG. 3: The phase diagram survival-extinction on (harmful 
mutation threshold T)- (mutation rate M) plane. 



of offsprings Bb shows some kind of scaling universality. 
Obviously, the larger mutational rate is the higher is the 
birth rate necessary to keep the population alive. On 
the other hand, the population size N(i) is restricted to 
the environmental capacity N max of the niche where this 
population lives: too large fertility of individuals may 
easily lead to species extinction due to overpopulation 
pfij . The overpopulation effect may be avoided when 
Verhulst factor acts only on the newly born babies 0] 
(see Fig. Etc)). 

In Fig. [3]thc phase diagram survival-extinction for var- 
ious values mutation rates M and the critical mutations 
threshold T is presented. The mutational pressure shifts 
damaged genes (positions of "1" in genome) in that way 
that only T — 1 mutations were located in the first R po- 
sitions of genome. This ensures chance of reproduction 
for individuals The population is stable for T > R 
independent of mutation rate M (see Ref. j2j|). For a 
lower threshold T < R the critical mutation rate M exists 
above which the population dies out due to mutational 
meltdown. 

The phase diagram on the plane of minimal reproduc- 
tion age R versus harmful mutation threshold T seems 
to show trivial behavior: when you allow individuals to 
reproduce immediately after they were born — the num- 
ber of individuals explode exponentially. When the pop- 
ulation size crosses the border of environmental capac- 
ity N max the Verhulst's factor will remove all individuals 
due to overpopulation. On the other hand when you shift 
the minimal reproduction age R almost to infinity (which 
here means N^u) the extinction will occur (see Fig. [3J 
as individuals will have no chance for reproduction. And 
again: the overpopulation effect may be avoided when 
one applies the Verhulst factor only on the infants. 

It seems that the transitions remain unchanged when 
the simulation time or the population size are enlarged. 
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(a) N max =10 6 , N =3 10 5 , B=4, M=1, E=32, m=b=1. 
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(b) N max =10 6 , N =3 10 5 , B=4, M=1, E=32, m=b=1. 
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FIG. 4: The phase diagram survival-extinction on (harmful 
mutation threshold T)-(minimal reproduction age R) plane 
for (a) Verhulst's factor acting on all individuals and (b) only 
on newly born babies. In the latter case the overpopulation 
is avoided. 



IV. THE CONCLUSIONS 

The subject of extinction is still attractive for scientists 
(see Refs. |2(| [2?j for the most recent papers in the 
field.) 

The influence of genetic factors on population extinc- 
tion was numerically investigated also earlier Q for the 
Eigen model Q of evolution, where the population docs 
not possess age structure. The Penna model with its set 
of control parameters is very flexible and thus almost any 
of experimentally observed effect for real population may 



be successfully reproduced . 

Some trivial assumptions of the values of the Penna 
model parameters may lead to species extinction. For 
example, when one assumes individuals infertility (B = 
b = 0) or one forgets to give individuals a chance for 
reproduction (E < 1 or R > Ntu) the initial population, 
even when huge, will die out. 

In this paper we show how the populations may be 
fragile when the main factor responsible for removing in- 
dividuals is either too rapid growth of the population 
(which results in exhausting resources necessary for life) 
or genetic mutational meltdown of the population. In 
both cases the phase transition extinction-survival is ob- 
served. When the environmental capacity is finite — and 
the Verhulst's factor guards maximal environment capac- 
ity border — even without any errors in reproduction (no 
mutations) the population may explode in one "season" 
much over allowed the limit and die out due to starva- 
tion. The effect of overpopulation may be removed when 
Vcrhulst factor kills only babies. On the other hand, for 
fixed average reproduction rate Bb the critical mutation 
rate Mm may be found above which population cannot 
have stable and non-zero size. 

The latter means, that even without catastrophe in 
the cosmic scale the last examples of the endangered 
species — for example the giant panda (Ailuropoda 
melanoleuca) , the red panda (Ailurus fulgens, better 
know as Firefox), Hickman's potentilla (Potentilla hick- 
manii), the San Francisco garter snake (Thamnophis sir- 
talis tetrataenia) j2^| — may be finally removed from the 
Earth and meet the fate of phoenix, dragons, gargoyles 
and trolls. But simulated extinction of such small popu- 
lations was already investigated in details earlier [171 ]. 

Last but not least: extinction not always must mean 
bad thing as the extinction of one species us ually open 
the free way for evolution of another one lll| . The 
most famous mass extinction — which killed all dinosaurs 
65 million years ago — perhaps was not very polite for 
these huge lizards (particularly from their point of view), 
but prepared place for a race to which author and anony- 
mous referee(s) of this paper belong to. 
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